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Abstract 

00 

In the Standard Model with Minimal Supersymmetry, the Lagrangian contains 
complex parameters which lead to additional CP violation. We study CP violating 
asymmetries in the decays of the MSSM charged Higgs boson H^, induced by loop 
corrections with intermediate SUSY particles, and perform analytical and numerical 
analysis. The decay rate asymmetry can go up to 25% and the forward-backward 
, asymmetry can reach up to 10%. 
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The CP-symmetry is broken in nature and this feature makes it an important tool 
for testing fundamental theories. An understanding of its non-conservation in particle 
physics may provide essential information about the mechanism, causing the dominance of 



P-i 

D | the matter over antimatter in the observable universe. Generally, effects of CP-violation 
r^, ' originate in non vanishing complex phases in the Lagrangian of the theory. In particular, 
in the Minimal Supersymmetric Standard Model (MSSM), the Higgs mixing parameter /i 
in the superpotential, the U(l) and SU(2) gaugino mass parameters Mi and M2, and the 
trilinear couplings Af (corresponding to a fermion /) can have physical phases which can- 
^ \ not be rotated away without introducing phases in other couplings [1]. The experimental 
upper bounds on the electron and neutron electric dipole moments (EDMs) constrain the 
phase of /i, M < O(10~ 2 ) [2] for a typical SUSY mass scale of the order of a few hundred 
GeV, however the phases of the other parameters mentioned above are practically uncon- 
strained. The CP- violating effects that might arise from the trilinear couplings of the first 
generation A U ( i are relatively small, as they are proportional to m u ^. However, the trilin- 
ear couplings of the third generation A t ^ yT can lead to significant effects of CP-violation, 
especially in top quark physics [3]. 

In the following, we study CP violation in the decays of the charged Higgs bosons H ± 
within the MSSM. At tree level, there are three possible decay modes of H + into ordinary 
particles: H + — ► tb, H + — > vt + and H + — > W + h°, where h° is the lightest neutral Higgs 
boson. Loop corrections due to a Lagrangian with complex parameters lead to non zero 
decay rate asymmetry between the partial decay widths of H + and H~, 
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and that would be a clear signal of CP violation. We consider such decay rate asymmetries 
in MSSM with complex parameters for the quark decay mode H + — > tb - the asymmetry 
o~ti P E9) an d f° r the bosonic decay mode H + — > W + h° p5] - the asymmetry 5^f h0 . For 
the quark decay mode H + — > tb we go a step further by including the decay products of 
the top quark. This allows to examine CP-violating asymmetries due to the polarization 
of the top quark [6] . The top-quark decays before forming a bound state due to its large 



mass, so that the polarization can be measured by the angular distributions of its decay 
products. Moreover, the polarization is very sensitive to CP violation. The considered 
CP violating forward-backward (FB) and energy asymmetries are constructed by using 
angular or energy distributions of the decay particles, following the formalism of [7] . This 
talk is entirely based on the results of [5] and [6]. 

We study the following processes, related to the quark decay mode of H ± [6] (See 

Fig.nD 
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The CP-violation is induced by SUSY loop corrections in the if ± t6-vertices, whose 
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Figure 1: The Feynman graphs of the processes we study. 

effective amplitudes are given by 

M H+ = iuipt) [Y b + P R + Y+P L }u(- Pl ) (4) 

M H - = iu{ Pb )[Y t -p R + Y b -P L ]U{- P t) , (5) 

The one loop diagrams that contribute are shown on Fig. [2j The Yukawa couplings 
(i = t, b) are a sum of the tree- level couplings y t and y b , and the contributions from the 
loops - 5Y i ± (i = t,b): 

Y± = y l + 5Y l ± i = t,b. (6) 

The loop contributions 5Y^ have, in general, both CP-invariant and CP-violating parts: 

8Y± = 5Yt nv ± \ 5Y t CP . (7) 

In addition, both CP-invariant and CP-violating parts have real and imaginary (absorp- 
tive) parts. Only the real parts Re SY t c b p enter the CP- violating asymmetries. We examine 
the following CP- violating asymmetries related to the processes (j2J) and Q: 




Figure 2: Sources for CP violation in H + — ► tb decays at 1-loop level in the MSSM with 
complex couplings (i, j — 1, 2; k — 1, 4). 



• CP-violating decay rate asymmetry 8 , defined as 

where A^j is the total number of particles /(/) in ([2]) or ([3]) respectively. 

• CP-violating forward-backward (FB) asymmetry AA^, which we construct 
from the FB asymmetries ± of the processes (J2J) and ([3]) using the angular dis- 
tributions of the fermions / = b', ^ from the top quark decay in (j2J) and . This 
asymmetry is defind as 

AAf =A™-A™, (9) 
where A± B are the ordinary FB asymmetries of the processes 

-pF _ -nB 

a fb _ 1 f,± 1 f,± nm 

1 f,± "T" 1 /,± 

and we have 

r£=/ - — ^cos#, rf=/ J—^d<x*e, (ii) 

J acoso 1 7| acoso 1 

i.e. are the number of particles /(/) measured in the forward direction of the decaying 
t(t) quarks and ^f± are the number of /(/) measured in the backward direction of the 
decaying t(t) quarks. 

• CP-violating energy asymmetry A.R^ P can be defined analogously, using the 
energy distributions of the processes (jSJ) and (j3J), namely 



Ai£ 6 P = R b+ — , (12) 
where R± are given by 

R r ± (x > Xp) -T±(x< Xq) . , 

^' ± "r ± (x>x ) + r ± (x<x )' 1 J 

x is a dimensionless variable proportional to the energy, x = E b > / and xo is any fixed 
value in the energy interval. 

The analythical results for the asymmetries are rather simple: 

tf P = *f P =S. (14) 



AAf p = ^mrnW H ™? + ™l )2 V CP , (15) 

U$ p = ± ab {m* H -rrS)'P GP , (16) 
where r CP and V CP are two linear combinations of the two formfactors: 

T CP = [y t Re(5Y t CP ) + (y b Re(SY b op )] (jp tPb ) - m t m b [y t Re(SY b cp ) + y b Re(5Y t CP 



(17) 



vCP _ y t Re{5Y t CP ) - y b Re(6Y b CP ) 



(Vt + y b )(PtPb) - 2m t m b y t y b ' 
Here T CP is the CP-violating part of the relevant H ± partial decay rate and V CP is the 
CP-violating part of the polarization of the top-quark. The coefficient a b m determines 
the sensitivity of the b-quark (lepton) to the polarization of the top quark: 



m? — 2m?,, 
+ 2m w 

Further analysis on the results shows that one does not need to explore numerically all 
of these asymmetries, as some of them are connected. In this context we reduce the 
numerical analysis on the following argumentation: 

• All of the discussed asymmetries are linear combinations of the CP-violating form 
factors Re(5Y t CP ) and Re(5Y b CP ), through a dependance on V CP or T CP (TTTD. 
In order to measure these formfactors, one needs to measure: 1) the decay rate 
asymmetries 5 b f which are proportional to T CP and 2) the angular and / or energy 
asymmetries which are proportional to V CP for a given process. 

• As there is no CP violation in t — > bW, the decay rate asymmetries in if ± — > W ± bb' 
and — > bb'l^v are equal to the decay rate asymmetry in iJ ± — > tb. We denote 
it, following @], by 5 CP : 

rCP r+ — r~ 

= 5 f p = 6f p = *— . 6 CP = — - , (20) 



where T^ 1 is the partial decay rate of the process H ± — > tb 



• The angular and energy asymmetries - AAff and AR b ' p , are not independent 
either. The angular asymmetries for the 6-quarks and the leptons are related by 

AAf p = ^ AA c b p « 2.6 AA c b p . (21) 

Furthermore, there is a simple relation between the 6-quark energy and angular 
asymmetries: 

A R CP — ( m H+ + m t) 2 A a cp ( 00 \ 

ARb ~ Am% + m* AA ' (22) 

which implies that for m#+ > m t , AR b p is bigger than AA^ P . Thus, in general, 
AR b p is the biggest asymmetry of those defined through the polarization of the top 
quark for m H + > 490 GeV. 

In the following, we present numerical results only on the decay rate asymmetry 5 CP and 
on the energy asymmetry AR b p . 

In order not to vary too many parameters we fix a part of the parameter space by the 
choice: 

M 2 = 300 GeV, M 3 = 745 GeV, M D = M Q = M b = M E = M L = 350 GeV, 
fx = -700 GeV, | A t | = \A b \ = \A T \ = 700 GeV. (23) 

According to the experimental limits on the electric and neutron EDM's, we take </> M = 

or M = 7r/10. The remaining CP- violating phases we vary, are the phases of A t , A b and 

A T . In the numerical code we use running top and bottom Yukawa couplings, calculated 

at the scale Q = m H +. Fig. [3^ and Fig. [3b show the asymmetries 5 CP and AR b p as 

functions of m H + for <p At = ir/2, (pA b = and </> M = 0. (Our studies have shown that 

the most important CP- violating phase is (f) At . There is only a very weak dependence on 

4>A b and 4>a and therefore we take them zero.) For tan/3 = 5 the decay rate asymmetry 
s cp 

goes up to 20%, while AR b p reaches 8% for the same values of the parameters. The 
asymmetries strongly depend on tan (5 and they quickly decrease as tan /3 increases. For 
mu+ < + m- b , the asymmetries are very small. The main contributions to both 5 CP 
and AR b p come from the self-energy graph with stop-sbottom. The vertex graph with 
stop-sbottom-gluino also gives a relatively large contribution. The contributions of these 
two graphs in 5 CP and AR b p are shown in Fig. H^i and Fig. Hb. The contribution of the 
rest of the graphs is negligible. 

Further, we take a very small phase, (j)^ = vr/10, in order to fit with the experimental data. 
As can be seen in Fig. Eh and Fig. [5b, the asymmetries can increase up to 25% for 5 CP 
and 10% for AR b p , respectively. The discussed asymmetries 5 CP and AR b p shows a very 
strong dependence on the sign of fi. As noted above, our analysis is done for fi = —700 
(see ([231 ). however if fi changes sign, fi = 700, all asymmetries become extremely small. 

We now turn to the bosonic decay mode of the charged Higgs boson H + W + h°. We 
are interested again in the CP violation caused by SUSY loop corrections in the H ± Wh°- 
vertex and study the decay rate asymmetry between the charge conjugate processes H + — > 
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Figure 3: The asymmetries 5 CP and AR^ P as functions of m#+ for </% = 7r/2, 0,4;, 
M = 0. The red, blue, and green lines are for tan/5 = 5, 10, and 30. 
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Figure 4: The contribution of the tb self-energy (red line), tbg vertex contribution (blue 
line) and the sum of the other (green line) diagrams to 5 CP and AR^ P as functions of 

m,H+ for tan/3 = 5 and (f>A t = 7r/2, 4>A b = 4>^ = 0. 



and H~ -> W^/i : 

CP r(#+ -> iy + /i°) - r(F- -> 



The matrix elements of if + — > W^/i and — > W~h° are given by 
M H ± = tge x a {p w )p a h Y ± 1 



(24) 



(25) 



where s^iPw) is the polarization vector of W , K are the loop corrected couplings. The 
CP- violating asymmetry is determined by the real parts Re (5Y h 
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5lf p = Re(5Yr ) + zlm(5rr ) 
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Figure 5: The asymmetries 5 CP and ARf p as functions of m#+ for <pA t = vr/2, <pA b = 
and a non zero phase of /i, M = n/10. The red, blue and green lines are for tan (3 = 5, 10, 
and 30. 

where the sum is over the loops with CP- violation shown on Fig. [6j We assume that the 
squarks are heavy and that the decay H + — > tb is not allowed kinematically. Thus, we 
are left with the loop corrections with sleptons, charginos and neutralinos. In the com- 
monly discussed models of SUSY breaking, the squarks are much heavier than sleptons, 
charginos and neutralinos. According to the numerical exploration of the — > W^h 
decay we have the following reasoning. 

First, increasing m#+, the mass m^o is saturated, approaching its maximum value, m^ x ~ 
130 GeV. There is an experimental lower bound for m^o, m h o > 96 GeV [8]. Thus, re- 
specting both the experimental and theoretical bounds, we consider m^o in the range 96 
< m h o < 130 GeV. Our analysis shows a very weak dependence on m h o and the results 
here are presented for m^o = 125 GeV. 

The second consequence concerns the H^Wh coupling which determines the BR (H + — > 
W + h°). It falls down quickly when increasing and, depending on tan /3, we can enter 
the so called decoupling limit, cos 2 (/3 — a) — > 0, where the BR (H + W + h°) almost 
vanishes. In order to keep the value of BR (H + — > W + h°) at the level of a few percents, 
we keep mn+ and tan/3 relatively small: 200 < < 600 GeV and 3 < tan/? < 9. 
In order not to vary too many parameters, we again fix part of the SUSY parameter space: 

M 2 = 250 GeV, M E = M L -5 GeV, M L = 120 GeV, \A T \ = 500 GeV, = 150 GeV. 

(27) 

In order 5£F to be nonzero, we need both new decay channels opened and CP-violating 
phases present. In accordance with this we have three cases: 1) When only the decay 
channels H + — > vf^ are open (Fig. [6k.). Then the phase <p T is responsible for CP viola- 
tion. 2) When the decay channels H + — > xtxt are °P en only (Fig. [6b). In this case CP 
violation is due to the phase 4> 1 , and 3) When both H + — > z>f+ and H + — > xtXk decay 
channels are kinematically allowed ( all diagrams of Fig. [6]). In this case the two phases 
r and 0i contribute. 

Examples of the asymmetry as function of for cases 1) and 2) are shown on Fig. [T^i 




Figure 7: Left: 5wh° as a function of a) for different values of tan/3, Ml = 120 GeV; 
solid lines are for <p T = —it/2, <pi = 0; dashed lines are for <p T = 0, <p\ = —ir/2. b) for 
different values of Ml, at tan/3 = 9, 4> T = — ir/2, 0i = 0. Right: S Wh o at tan/3 = 9 
versus the CP violating phases T and X for different values of M [ M = 0, 7r/4, 7r/2]. 
The solid lines are for r = [— 7r, it] while 0! = 0; the dashed lines are for <pi = [—it, it] 
while (f> T = 0, a) for m H + = 237 GeV (= rriy + m~+) and b) for m H + = 387 GeV 

(= m _+ +m -o). 



for different values of tan (5. It is clearly seen that in both cases the asymmetry strongly 
increases with tan/3. In both cases, 1) and 2), the asymmetry reaches up to 10 -2 . The 
dependence on Ml for case 1) is seen on Fig. [7b. Case 3), when all relevant SUSY par- 
ticles can be light, is described by the algebraic sum of the two graphs at a given tan/? 



and we don't present it separately. In all these cases the asymmetry does not exceed a 
few per cents. 



The effect of a non-zero phase M is seen on Fig. [71 In all cases a CP violating phase 
of \x does not change the form of the curves but rather shifts the positions of the maxima 
and, in general, increases the absolute value of the asymmetry. 

Let us summarize. We have calculated different asymmetries concerning the quark 
and boson decay mode of the charged Higgs boson H ± , caused by CP- violating phases in 
the MSSM Lagrangian. The quark decay mode is dominant for large and the most 
important phase is the phase of A t . In this case the decay rate asymmetry can reach up 
to 25%, and the forward-backward asymmetry can go up to 10%. These asymmetries are 
rather large and in principle measureable at LHC, but because of the large background of 
the process, require higher luminosity (e.g. SLHC). The boson decay mode is important 
for small m#+ and tan j3 and sensitive to the phases of Mi and the phase of A T . The decay 
rate asymmetry is, typically, of the order of 10~ 2 -r- 10~ 3 . Concerning its measurability 
more promising are the next generation linear e + e~ colliders. 
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